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Abstract 

In this paper, an n-user Gaussian interference channel, where the power of the 
transmitters are subject to some upper-bounds is studied. We obtain a closed-form 
expression for the rate region of such a channel based on the Perron-Frobenius theorem. 
While the boundary of the rate region for the case of unconstrained power is a well- 
established result, this is the first result for the case of constrained power. We extend 
this result to the time-varying channels and obtain a closed-form solution for the rate 
region of such channels. 

I. Introduction 

Channel sharing is known as an efficient scheme to increase the spectral efficiency 
of the wireless systems. While such a scheme increases the capacity and the coverage 
area of systems, it suffers from the interference among the concurrent links (co-channel 
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interference). Consequently, the signal-to-interference-plus-noise-ratio (SINR) of the 
links are upper-bounded, even if there is no constraint on the transmit powers. 

There have been some efforts to evaluate the maximum achievable SINR in the 
interference channels. In [1], the maximum achievable SINR of a system with no 
constraint on the power is expressed in terms of the Perron- Frobenius eigenvalue of a 
non-negative matrix and this result is utilized to develop an SINR-balancing scheme 
for satellite networks. This formulation for the maximum achievable SINR is deployed 
in many other wireless communication applications such as [2]-[5] afterwards. 

Recently, the rate region of interference channels and its properties has been 
investigated in the literature. In [6], it is shown that the capacity region when the 
power is unbounded is convex. The capacity region in [6] is defined as the set of feasible 
processing gains while for a constant bandwidth, the processing gain is inversely 
proportional to the rate. In [7], some topological properties of the capacity region 
(with the aforementioned definition) of CDMA systems are investigated for the cases 
when there are constraints on the power of individual users and when there is no 
constraint on the power. The authors in [7] show that the boundary of the capacity 
region with one user's power fixed and the rest unbounded is a shift of the boundary 
of some capacity region with modified parameters, but unlimited power. However, 
this result is not in a closed form and can not be extended for the other forms of 
power constraints. 

It is shown that the feasible SINR region is not convex, in general [8]-[10]. In 
[11], it is shown that in the case of unlimited power, the feasible SINR region is log- 
convex. The authors in [6] also consider a CDMA system without power constraints, 
and show that the feasible inverse-SINR region is a convex set. In [8], it is proved 
that the feasible quality of service (QoS) region is a convex set, if the SINR is a 
log-convex function of the corresponding QoS parameter. Reference [12] shows that 
under a total power constraint, the infeasible SINR region is not convex. 

In this paper, we extend the result on the maximum achievable SINR in [1] 
to the systems with certain constraints on the power of transmitters. This result 
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which is based on Perron-Frobenius theorem, yields a closed-form solution for the 
rate region of the systems with constrains on the power. The extendable structure 
of constraints enables us to use the proposed derivation for the maximum achievable 
SINR in many practical systems. This result is extended to a time-varying system, 
where the channel gain is selected from a limited-cardinality set, and the average 
power of users are subject to some upper-bounds. 

Notation: All boldface letters indicate column vectors (lower case) or matrices 
(upper case). Xij and Xj represent the entry and column i of the matrix X, 

respectively. A matrix X nxm is called non-negative if > 0, Wi,j, and denoted by 
X > 0. Also, we have 



where X, Y and are non-negative matrices of compatible dimensions [13]. det(X), 
Tr(X), X', and |X| denote the determinant, the trace, the transpose, and the norm of 
the matrix X, respectively. I is an identity matrix with compatible size. <S> represents 
the Kronecker product operator, diag(x) is a diagonal matrix whose main diagonal is 
x. We define the reciprocal of polynomial q(x) of degree m as x m q(^). ip(X,y,S) is 
a matrix defined as a function of three parameters, which are respectively a matrix, 
a vector and a set of indices, 



An interference channel, including n links (users), is represented by the gain 
matrix G = [gij] n xn where g^ is the attenuation of the power from transmitter j to 
receiver i. This attenuation can be the result of fading, shadowing, or the processing 
gain of the CDMA system. A white Gaussian noise with zero mean and variance of is 
added to each signal at the receiver i terminal. In many applications, the QoS of the 
system is measured by an increasing function of SINR. In an interference channel, 



X > Y 



X- Y > 0, 




II. Problem Formulation 
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SINR of each user, denoted by %, is 



7»= ^? > Vi G {!,..., n}, 



where pi is the power of transmitter i. In addition, in practice, the power vector p is 
subject to a set of constraints. The main goal is to find the maximum SINR which 
can be obtained by all users in the presence of such constraints. To this end, we solve 
the following optimization problem 

max 7 (1) 
s.t. ji > [id (2) 
p > (3) 
J2Pi<Pn, (4) 

where Q C {l,...,n} with k elements and n is a given vector with fa > and 
= 1. As we will see, the solution can be easily extended for the case of multiple 
power constraints of the form ^^p? < Pq for different f2 C {l,...n}. fi provides 

the flexibility of satisfying different rate services for different users. According to 
Fig. [TJ the solution of ([T|) yields the maximum achievable SINR in the direction of 
vector \x. Although the numerical solution of this problem is already obtained through 
geometric programming [14], [15], we propose a different approach which leads to a 
closed- form result. 

By defining the normalized gain matrix A as 



A — [O-jj'Jnxni &ij — \ ^ i 



9ij ■ , ■ 

* T 3 



the constraint ([2]) is rewritten as 

Pi 



i = 3 



> 7 , Vie{l,...,n}, (5) 



Fig. 1. The boundary of SINR Region for an interference channel with 2 users 



where 

Vi = - JlJl , V = [Vi]nxl- (6) 

9a 

Since we are interested in maximizing the minimum SINR, if SINR of one user is 
more than that of the others, it can reduce its power to other users' advantage, and 
finally the minimum SINR is improved. Therefore, equality holds in (JHJ) as 

P = 7, V*G{l,...,n}. 

3=1 

After reformulating the problem in a matrix form we will have 

(-I-diag(/^)A)p = 77. (7) 

7 

The objective is to find the maximum 7 while the system of linear equations in (jTj) 
yields a power satisfying the constraints on the power vector fl3J), (jlj). 

When there is no constraint on the power vector (rather than trivial constraint 
of p > 0), the maximum achievable SINR, 7*, is characterized based on the Perron- 
Frobenius theorem as 

1 

7 



A*(diag(/x)A) ' 
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where A* is the Perron- Frobenius eigenvalue of the associated matrix [13]. This result 
was deployed in the communication systems for SINR-balancing (fii = /i 2 = . . . = fi n ) 
in [1] for the first time. 

We find the maximum achievable SINR, considering certain upper-bounds on 
the power of transmitters in the following sections. 

III. SINR Region Characterization 

We define F as 

F = I - 7 diag(^)A. (9) 

Then, the system of linear equations in ([7]) is reformulated as 

Fp = 777, (10) 

where rj is defined in ([6]). According to the Cramer's rule, the solution to ffTUl) is 
obtained by 

det(HW) 

Pi - 



det(F) ' 

where 



H« = [h«] nxn , h? = \ ™ j ~ l . (II) 

f i 3 + 1 



Defining h.w( 7 ) = det(H^) and f( 7 ) = det(F), we have 

h(%) 

Therefore, the constraint in (HI) can be written as 



£h W (7) 



f(7) 
Defining 

un(7)=^f(7)-E h(i) W' 



< Pn- (12) 
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the inequality f|T2|) is equivalent to 

Uo(t) 



>0. (13) 



f(7) 

We desire to find the largest possible interval where both the numerator and the 
denominator have the same sign. It can be shown that this interval is connected and 
adjacent to zero. Apparently, un(0) > 0, and f(0) > 0. Consequently, 

3 e > : f(e) > and u n (e) > 0. 

Therefore, both the numerator and the denominator are positive in the positive 
neighborhood of zero. For satisfying ([TBI , we have to find the smallest positive real 
simple root of the numerator and the denominator, r(un) and r(f), and take the 
minimum of the two as 

7 = min{r(f),r(un)}. (14) 

For the sake of simplicity, without loss of generality, we assume that O = {1, . . . , k}, 
k < n, i.e., the first k users are subject to the total power constraint. For the 
numerator we have 

k 

un(j)=p n det (F)-^det (H») 

i=l 
k 

= ^(det(F)-^det(H«)), (15) 

i=l 

where flw is defined as 

— ' = i 
f i 3 ^ i 

Lemma 1 If square matrices X and Y differ only in column i, i.e., 

*j ^ Yj J = i 
x i = Yj J + i 
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then 

det (X) + det(Y) = det (^(X, y h {«})) 
= det (</,(Y >Xi , {«}))• 

Equation (|T5j) is rewritten as 

fe 

un(7) = Pn(det (F) - det(H«) - ^det (H«)). (16) 

Since F and IT 1 ** are the same except for the first column, using Lemma [I] , we will 
have 

det (F) - det(H (1) ) = det (^(F, {l})). (17) 

Pn 

On the the other hand, using the fact that addition or substraction of columns does 
not change the value of the determinant, we will have 

det (H«) = det (V(H (i) , -hf ,{1, 1})). (18) 

Then, using (fl7l) and (fl8l) and regarding h- = we can rewrite ([TBI) as 



un(7)=Pb(det(^(F,-^,{l})) (19) 
v Pn 

-Vdet^(H®,-^,{l,...,i-l}))). 

i=2 P ^ 

Since F and are the same except for the column z, we can easily see that 

■^(F, — {1, . . . , i — 1}) and ^(H", — • = — , {1, . . . , % — 1}) are the same except for 

Pn 

the i column. Therefore, 

det (^(F,-^,{l,...,z-l})) 
Pn 

-det (^(H«,-^,{l,...,i-l})) 
Pn 

= det (,/,(F,-^,{l, ...,*})). 
Pn 

Applying this result to (fl9|) successively yields the following lemma. 
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Lemma 2 

110(7) =Pn det ty,(F,-^,fi)). 

Pn 

We utilize the result in Lemma [2] to find the smallest positive simple root of uq using 
Perron- Frobenius theorem. This theorem states some properties about the eigenvalues 
of a primitive matrix. A square non-negative matrix X is said to be primitive if there 
exists a positive integer k such that X fc > [13]. 

Theorem 1 [13] (The Perron- Frobenius Theorem for primitive matrices) Suppose 
X is an m x m non-negative primitive matrix. Then there exists an eigenvalue A*(X) 
(Perron-Frobenius eigenvalue or PF- eigenvalue) such that 

(i) A*(X) > and it is real. 

(ii) there is a positive vector v such that Xv = A*(X)v. 

(iii) A*(X) > |A(X)| for any eigenvalue A(X) ^ A*(X). 

(iv) //X > Y > ; then A*(X) > |A(Y)| for any eigenvalue ofY. 

(v) A*(X) is a simple root of the characteristic polynomial o/X. 

Lemma 3 The smallest positive root of uq(j) is 

r(im) = — r) r- 

A*(^(diag(//)A != ^)) 
Pn 

Proof: 



Ml)=Pn^ (V(F,-^,fi)) 

Pn 



Pn det (if)(I - 7diag(/x) A, -— , fi)) 

Pn 

p n7 n det (t/>(-I - diag(^) A, -I, O)) 

7 Pn 

p c7 n det(-I-^(diag(^)A,^-,fi)). 

7 Pn 
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Consequently, — is the reciprocal of the characteristic polynomial of the matrix 

PnT 

■0(diag(//)A, — ,ft). Therefore, the roots of this polynomial are equal to the inverse 
Pa 

of the eigenvalues of ■0(diag(^)A, — ,Q). On the other hand, according to Theorem 

[H since ^(diag(/j)A, — ,Q) is a primitive matrix, the PF-eigenvalue of this matrix 

Pn 

is real and positive and has the largest norm among all eigenvalues. Also it is the 
simple root of the characteristic polynomial of the aforementioned matrix. Therefore, 
the inverse of this eigenvalue gives the smallest positive simple root of un(7) and the 
claim is proved. ■ 
For the denominator using (Q, we have 

f( 7 ) = det (F) = det (I - 7 diag(/j)A) 

= 7 " det (-1 - diag(/x) A) . (20) 

7 

Therefore, f( 7 ) is the reciprocal of the characteristic polynomial of diag(/x)A. On 
the other hand, according to Theorem (TJ the PF-eigenvalue of diag(/j)A, is real 
and positive. It also has the largest magnitude (norm) among the eigenvalues of the 
matrix and it is the simple root of the characteristic polynomial of the associated 
matrix. Therefore, A*(diag(/j)A) is the inverse of the smallest positive simple root 
of diag(/j)A. Thus, 

r(t) = —7 r- (21) 

v ; A*(diag(/j)A) v ; 

On the other hand, according to (jSJ), r(f) is also the maximum achievable SINR 

for the system with unbounded powers satisfying constraint ([3]). Consequently, using 

(I14p . (|2T|) and LemmaQ, the maximum achievable SINR to satisfy all constraints 

on the power (constraints (J3]) and (j3J)) is 

7* = min {r(f), r(u^)} 

= mm{ — -, — - — r-, Yi }• 

A* (diag(/j) A) x * (V»(diag(/i) A, ±, {*})) 

Pi 
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Since ■0(diag(/i)A, — , flj > diag(/x)A and both are primitive, using Theorem [1] we 

Pi 

have 

A*(#Jiag(j*)A,?,{i})) > A*(diag(/x)A), 

Pi 

and consequently the maximum achievable 7 for a system with constraint on the total 
power of any subset of the users is achieved. This discussion leads to the following 
theorem. 

Theorem 2 The maximum achievable 7 in an interference channel with n links and 
gain matrix A ; where power vector is subject to the following constraints, 

P>0, 

^2 Pi < Pn 

ien 

is equal to 

• = 1 

A*(^(diag(/*)A,^-,fi))' 
Pn 

where Q C {1, . . . , 71} is an arbitrary subset of the users. 

When multiple constraints on power exist, it is obvious that the maximum 
achievable SINR is the minimum of the maximum achievable SINR when each of 
the constraints is applied separately, i.e., 

7 * = min7*, (22) 

where 7* is the maximum achievable SINR for the constraint i on power. The following 
corollary yields the maximum achievable SINR when the power of individual users 
and the total power are constrained. 
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Corollary 1 The maximum achievable 7 in §T§, where power vector is subject to the 
following constraints, 

P>0, 
P < P, 

n 

^2pi <p t 

i=i 

is equal to 7* = 

min{ — -, (23) 

A*(^(diag(/x)A,^,{l,...,n})) 

1 

A*(#liag(/x)A,^,{l}))' 
1 

A*(^(diag(^)A,^-,{2}))' 
P2 

I } 

A*(^(diag( M )A,^-,{n})) 

Pn 

The boundary of the SINR region in any direction can be obtained by choosing 
/x, accordingly. Due to the explicit relationship between the SINR and the rate in 
Gaussian channels, obtaining the SINR region in these channels amounts to the rate 
region characterization. As an example, Fig. [2] and [3j respectively, depict the rate 
region and SINR region of a system with the gain matrix G as 

0.6791 0.0999 

G = 

0.0411 0.6864 

while the power of individual users and the total power are upper-bounded as p 1 = 
0.8, p 2 = 1, p t = 1.4, and o\ = a\ = 10 _1 . 

The rate region is simply the intersection of all the rate regions resulted from 
applying each constraint separately. As shown in Fig. [2] and Fig. [HI the boundary 
of SINR and rate region, when there is no upper-bound on powers is always above 
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■ A 

■ B 

c 

• D 



Fig. 2. The rate region for a 2-user interference channel with the following constraints on the power, 
A: pi >0,p2> 0, B: pi + p 2 > Pt, Pi > 0, p 2 > C: < pi < pi, p2 > 0, D: < pi < p2, pi > 




Fig. 3. The rate region for a 2-user interference channel with the following constraints on the power, 
A: pi >0,p 2 > 0, B: pi + p 2 > pt, Pi > 0, p 2 > C: < Pl < pi, p 2 > 0, D: < p 2 < p 2 , Pi > 
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012345678 



Fig. 4. The rate region for a 2-user interference channel with the following constraints on the power, 
A: pi >0,p2> 0, B: pi + p 2 > pt, Pi > 0, p 2 > C: < pi < pi, p2 > 0, D: < pi < p2, pi > 




Fig. 5. The rate region for a 2-user interference channel with the following constraints on the power, 
A: pi >0,p 2 > 0, B: Pl + p 2 > pt, Pi > 0, p 2 > C: < pi < pi, p 2 > 0, D: < p 2 < p2, Pi > 
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other boundaries. It is because of the fact that the maximum achievable SINR for the 
unbounded-power system is the inverse of PF-eigenvalue of diag(/x) A; while the maxi- 
mum achievable SINR when the power is bounded, is the inverse of PF-eigenvalue of a 
matrix which is definitely greater than diag(/x)A. Therefore, based on Theorem [1] the 
unbounded SINR boundary would be above the bounded-power systems. Thus, this 
boundary doesn't have any direct role in forming the main boundary. An interesting 
observation is that if the p^s or p t are increased the boundaries of bounded-power 
systems tend to the unbounded-power system boundary; the extreme case is when 
the maximum power goes to infinity which means the power is unbounded, then the 
matrices whose inverse of PF-eigenvalue form the boundaries become equal and these 
boundaries touch each other. 

As another observation, the rate and SINR regions for a 2-user channel with 
weaker cross links are shown in Fig. H] and [H The gain matrix in this system is 
assumed to be 

' 2.0430 0.0359 
0.0134 1.3313 



G 



(24) 



while the power of individual users and the total power are upper-bounded as p 1 = 

1, p 2 = 1, % — 1-5, and a\ = a\ = 10 _1 . The extreme point of this situation is 

when the links have no interference on each other and therefore the maximum SINR 

for each user considering the individual constraints would be SINRi = — — ^ — for 

a 2 

each user. We can see in Fig. [4] and Fig. [5] that these boundaries are more straight 
than the ones in Fig. [2] and Fig. [3] which confirms our conjecture. 

IV. Time- Varying Channel 

So far, we have assumed that the channel gains are fixed with time. However, in 
practice, channel gains vary with time due to the users' movement or changing the 
environment conditions. 

In this section, we consider an interference channel with n co-channel links 
whose channel gain matrix is randomly selected from a finite set {Gi, . . . , G/} with 
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probability pi, . . . , pi, respectively. The matrix Aj denotes the normalized gain matrix 
in the state i, i G {1,...,/}. The objective is to find the maximum 7 which is 
achievable by all users in all channel states, while the average power of the users are 
constrained, i.e., 

max 7 

S.t. 7j+(i-i) n > fijj, Vj 6 O, i e {1, ■ ■ ■ , 1} 

Pj+(»-i)n > 0, Vj e n,i e {i,...,0 (25) 

where 7j + (j_i)„ and Pj + (j_i) n are the SINR and the power of transmitter j respectively, 
when the channel gain matrix is Gj. We define an expanded system including In users 
with block diagonal matrices G and A as the channel gain matrix and the normalized 
gain matrix, respectively. In the matrices G and A, the i th matrix on the diagonal 
is Gj and Aj, respectively. It is clear that block diagonal format of these matrices 
indicate that there is no interference between the links associated with different states. 
Like the previous discussions the requirements on these links form a system of linear 
equations with the following formulation in a matrix form, 

(-hnxin ~ diag(li X i <g> /x)A)p = T], (27) 

7 

where 

2 

Vj+(i~i)n = ,j ett,i e {l,. . . ,1}. 

9j+(i-l)n,j+(i-l)n 

According to we define F as 

F = -linxin - diag(l ix i <g> /x)A. 

7 

Then, we have 

Fp = 777. 
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Using Cramer's rule, we will have 



Pj+(i-l)n 



det(F) 

where H <J+< - J_1 ) n ) according to (TTTT) is the matrix F whose column j + (i — l)n is 

substituted by 777. The average of the total power of the users in Q is equivalent to 

1 

v det(H^+^- 1 H) 
i=i j&n v y 

^E^E det ( H(i+( '" 1)n) )- ( 28 ) 



det(F / 



Based on ([28]) . we define 



1 

)n)- 



uo( T ) = p n det(F) - X> £ det(H^+( 1 - 1 

i=l j'GO 

and 

f( 7 ) = det(F). 

Therefore, the constraint in (1261) is equivalent to 

un(7) > n 
f(7) " ' 

Like before, it is easy to show that the maximum achievable SINR satisfying con- 
straints (J25I) and flUD is 

7* = min{r(f),r(u n )}. (29) 

To simplify u^(7), we have 

1 

un( 7 ) = det(F) - £ p,^ detCH^*" 1 )")) 

i=i jen 

= Pn(det(F) - £ £ det^^)), 

i=l jefi 
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where H( J+ (* ^ is H^ J ' + ^ ^ whose column j + (i — l)n is multiplied by Using 

Pn 

the same procedure as before, we obtain 

u n (7) =Pndet (F-D), 

where 

i 

D = ^ ^(0 faxfa , {j + (i - l)n : J G fi}) 

i=i p « 

It is easy to see that 

r (un) = 1 • 

A*(diag(l* x i <8> n)A + V ^(Oi„xin, {j + (* - l)n : j G fi})) 

<=i p ^ 

and 

r(f) ' 



A*(diag(l ixl ®M)A)' 
Therefore, using Theorem [1] and equation (|29]1 . we will have the following theorem. 

Theorem 3 The maximum achievable 7 in a time-varying interference channel with 
n links and probability vector p ixl , with the following constraints on power, 

Pj,i > 0,Vj G tt,i G {1,...,/}, 

is equal to 

1 



7 



A*(diag(l ixl ® fi)A + V V(O tax «n, {j + l)n : j G 0})) 

i=i p « 

Apparently, if there are multiple constraints on the power, the maximum achievable 
SINR 7* is computed by 

7* = minT*, 
i 

where 7* is the maximum achievable SINR obtained by Theorem [3] while only the 
constraint % is considered for the system. 
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V. Conclusion 

In this paper, we have obtained a closed-form solution for the maximum achiev- 
able SINR in an interference channel, utilizing the Perron- Frobenious theorem, when 
there is a total power constraint on any subset of the users. This result leads to 
characterizing the boundary of the rate region with multiple constraints on the power. 
In addition, we considered a time-varying interference channel where the average of 
total power of an arbitrary subset of the transmitters is subject to an upper-bound. A 
closed-form expression for the rate-region of such a channel is obtained and extended 
to the systems with multiple power constraints. 
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